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We study the effective quark mass induced by the finite separation of the domain walls in the domain-wall 
formulation of chiral fermion as the function of the size of the fifth dimension (L s ), the gauge coupling (/3) and 
the physical volume (V). We measure the mass by calculating the small eigenvalues of the hermitian domain- wall 
Dirac operator (-Hdwf(tio = 1-8)) in the topologically-nontrivial quenched 517(3) gauge configurations. We find 
that the induced quark mass is nearly independent of the physical volume, decays exponentially as a function of 
L s , and has a strong dependence on the size of quantum fluctuations controlled by /3. The effect of the choice of 
the lattice gluon action is also studied. 



Simulating massless or near-massless fermions 
on a lattice is a serious challenge in numerical 
quantum field theory. The origin of the diffi- 
culty can be traced to the well-known no-go the- 
orem first shown by Nielsen and Ninomiya, which 
states that one can not write down a local, hermi- 
tian, and chirally-symmetric lattice fermion ac- 
tion without the fermion doubling problem [Q. 
Hence to have chirally-symmetric fermions on a 
lattice, one must use nonlocal actions in which 
the coupling between lattice sites do not identi- 
cally vanish even when the separation becomes 
large, as long as one insists on having a 4 dimen- 
sional action. 

One of the lattice chiral fermion formalisms 
that have been studied extensively in recent years 
is the domain- wall fermion, first formulated by 
D. Kaplan Q and later modified for realistic lat- 
tice simulation by Shamir ||. In the domain- 
wall construction, one introduces an extra fifth 
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dimension s with a finite extension. After dis- 
cretization, the fifth direction has L s number of 
lattice sites. If we put the same four-dimensional 
gauge configuration on every four-dimensional s 
slices, the five-dimensional massive theory admits 
a four-dimensional effective theory in which a left- 
handed chiral fermion lives near the s = slice 
and a right-handed one near s = L s — 1. In- 
tegrating out the heavy modes [Q, one obtains 
an effective four-dimensional chiral theory in the 
limit of L s — > oo. For finite L s , however, the two 
chiral modes can couple to produce an effective 
quark mass. Strong gauge field fluctuations can 
induce rather strong coupling, and hence rather 
large quark mass. This quark mass is expected to 
decrease exponentially as L s — > oo with possible 
power law corrections. 

In Ref. (5), we proposed to measure the in- 
duced effective quark mass by considering the 
eigenvalue of the hermitian domain- wall Dirac op- 
erator HuwF(m ) = R s "f5DuwF(m Q ), where R s 
denotes the reflection in the fifth direction. In the 
L s — > 00 limit, the lattice version of the Atiyah- 
Singer theorem |6) guarantees the existence of ex- 
act zero modes in the presence of a instanton 
background. For finite L s , the lowest eigenval- 
ues of the Dirac operator are not zero. We take 
the average of these would-be-zero eigenvalues as 
the effective quark mass. The result is qualita- 
tively consistent with those obtained from other 
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Crossing points vs. Effective masses (m - 1.8) 

L =8,12,16, m ta = 21.2(2.8),12.0(2.4),8.1(2.1)xl0 -3 , f= 5.85, 8 4 lattice 
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Crossing point vs. Effective mass (m Q =1.8) 

L =16,20,24, m =17.8(1. 0),13.9(0.9),11.4(0.8)xl0~ 3 . p=5.7, 8 4 lattice 
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Figure 1. Effective quark mass induced by do- 
main walls with domain wall height mo = 1.8 for 
the Monte Carlo configurations at (3 = 5.85, 8 4 
lattice. L s is the number of lattice sites in the 
fifth direction. 



Figure 2. Effective quark mass induced by do- 
main walls with domain wall height mo = 1.8 for 
the Monte Carlo configurations at (3 — 5.7, 8 4 lat- 
tice. L s is the number of lattice sites in the fifth 
direction. 



methods j7|-§. Here we report a more systematic 
study of the effective quark masses to understand 
the effects of different lattice sizes, the coupling 
constant /?, and the form of gluon actions. 

To study the effect of the gauge coupling, we 
have generated 100 5/7(3) lattice gauge configu- 
rations on a lattice size 8 4 , 50 each at /3 = 5.85 
and 5.7. We measure the spectral flow of the 
hermitian Wilson-Dirac operator to calculate the 
topological index of the gauge configurations, and 
calculate the eigenvalues of 7?dwf(wo = 1.8) cor- 
responding to the nontrivial topology of the gauge 
configurations. For larger lattice spacing, quan- 
tum fluctuations are stronger, and some of these 
fluctuations can be misidentified as small size in- 
stantons. It turns out that they can induce strong 
couplings between the left-and right-handed chi- 
ral modes and are detrimental to the existence of 
the low-energy effective theory. Indeed, for the 
same L s , the effective masses are much larger at 
the smaller /3's than those at (3 — 6.0. For exam- 
ple, with L s = 16, m off is 0.00080(13) at (3 = 6.0, 
0.008(2) at (3 = 5.85, and 0.018(1) at j3 = 5.7. 

Moreover, for the same lattice size, physical 
volume is larger at smaller f3, and hence can house 
more fermionic zero modes. As shown in Fig. [j] 



for (3 = 5.85 and Fig. | for (3 = 5.70, the total 
number of fermionic zero modes in the 50 config- 
urations is now 32 and 96, respectively. 

For (3 = 5.7, there are several crossings very 
close to m = 1.8. This is in contrast to what 
has been observed at (3 = 6.0, where the cross- 
ings occur mostly around m — 1.0 PJlC|]. The 
equal spacing (in a log plot) between m e s at dif- 
ferent L s is a clear signal for the exponential de- 
cay. However, there is a significant variation in 
the rate among all the crossings, as evident in the 
figures. (Note: The recent works of the CP-PACS 
and RBC || collaborations show the signal 
for varying rate of exponential decay and/or non- 
vanishing effective mass in the L s — > oo limit. 
This behavior is only seen at a much larger L s 
than those studied here. However, it is interesting 
to note that averaging eigenvalues with varying 
exponential rate can easily reproduce the large 
L s behavior of the effective mass observed in the 
aforementioned reference.) 

To see the volume dependence at a fixed mo 
and (3, we also measure the effective mass on a set 
of 50 configurations on an 8 3 x 16 lattice at (3 = 
5.85. The total number of fermionic zero modes 
is now 64, doubling that on the 8 4 lattice. The 
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Crossing points vs. Effective masses (nip= 1.8) 

L =8,12,16, m =15.9(0.6), 8.7(0.4),6.0(0.4)xl0~ 3 , B= 7.9, 8 x16 lattice 
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Figure 3. Effective quark mass induced by do- 
main walls with domain wall height tuq — 1 .8 for 
the 50 Monte Carlo configurations generated by 
the 1-loop, tadpole improved gauge action |ll[] on 
a (3 = 7.9, 8 3 x 16 lattice. L s is the number of lat- 
tice sites in the fifth direction. 



Figure 4. Average coefficient of the exponential 
decreases as a function of p(0;mo). The extrap- 
olation of the fit to the continuum limit for the 
Wilson gauge backgrounds is shown on the left. 
p(0; mo) for the Wilson and improved gauge ac- 
tion are from Ref. |l(J for 8 3 x 16 lattices. 



average effective mass turns out to be essentially 
the same as that on the smaller lattice. 

We also measure the effective mass on a set 
of 200 configurations generated from the 1-loop 
Liischer-Weisz gauge action ]ll| with the tadpole 
improvement. (Crossings from only 50 lattices 
are shown in Fig. ^|) Similar studies using var- 
ious improved gauge actions, including an RG 
improved action [jl2| are reported in 0,|l3|- The 
gauge coupling {(3 = 7.9) corresponds to a spac- 
ing of ~ 0.16 fm, similar to (3 — 5.7 of the Wilson 
action. The spectral flow and the domain-wall 
eigenvalues are studied with the same Wilson- 
Dirac operator. As shown in Fig. ^, the num- 
ber of instantons as well as the distribution of 
the crossings differs significantly from the Wilson 
action at similar lattice spacing. Because of the 
decrease of small-scale quantum fluctuations, the 
probability of the crossing at mo > 1.2 is heavily 
suppressed, and the density of small eigenvalues 
of Hw(m ) is also much smaller (<~ 2.7 x 10~ 4 , 
compared to ~ 1.8 x 10~ 3 for (3 = 5.7 Wilson 
action). Therefore, the effective mass decreases 
faster as a function of L s . The average effective 



mass at L s = 16 is ~ 6 x 10~ 3 , compared to 
18 x 10~ 3 from the Wilson action at (3 — 5.7. 

This point is also reflected in the dependence of 
the exponential decay rate on the density of the 
zero eigenvalues of Hyj{m ) (p(0;mo)). We note 
that since the gauge fields are replicated along 
the fifth dimensional slices, the relevant length 
scale in the fifth dimension is the inverse of the 
rate of exponential decay (a) and by simple en- 
gineering dimensions is given qualitatively by the 
density of zero eigenvalues of p(0; mo) 1 / 3 . In Fig. 
[|, we have plotted 1/a as a function of p(Q;mo). 
The rate a for each coupling is calculated by fit- 
ting m g{L s ) = moexp[— aL s ] at different L s . 
The statistical errors are estimated by doing cor- 
related fits to single-eliminated jackknife blocks. 
The inverse decay rates from all the configura- 
tions studied show an approximate linear scaling 
as p^mo) 1 / 3 . This suggests that the density of 
small eigenvalues is indeed the dominating fac- 
tor for the exponential decay rate of domain-wall 
fcrmions. 

Effective masses thus obtained for the differ- 
ent gauge coupling and L s are plotted in Fig. ||. 
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Average effective masses 




G-0 5.7, 8 x32, GMOR 
B-B 5.7, 16 3 x32, GMOR 
<M>5.7,S 4 , topology 
&r-& 5.85, S 4 , topology 
<^ 5.85, 8 3 X16, topology 
V~V 6.0, 8 , topology 
X— X 6.0, 16 3 x32, 
Jf^-K LW, 8 3 X|6, topology 
H h 6.0, 16 3 x32, WT 



sense expected because the coupling of the chiral 
modes between the opposite walls has little to do 
with the size of the four-dimensional slice. On the 
other hand, a strong dependence on f3 is observed. 
In particular, the effective mass is much larger at 
P = 5.85 or 5.7 than that at j3 = 6.0. For the 
improved gauge action, the spurious fluctuations 
are reduced significantly and the L s needed to ob- 
tain a good chiral symmetry is reduced. Since the 
additional computation needed for the improved 
action is negligible, using improved gluon actions 
may enable us to simulate domain-wall fermions 
with larger lattice spacing. 



Figure 5. Average effective masses from vari- 
ous observables as a function of L s . Effective 
masses from (3 = 5.7, GMOR relation are from 
Ref. p4| . Data for (3 — 6.0, and the ax- 
ial Ward-Takahashi(WT) identity are from Ref. 
H and 0, respectively. LW denotes the 1-loop, 
tadpole improved gauge action from Ref. pTf . 



The data from Ref. are also included for 

comparison. For a given gauge coupling, differ- 
ent volumes and methods of measurements have 
little effect on the size of the effective mass as 
well as the rate of exponential decay. However, 
the change of gauge action affects the effective 
mass significantly. It is quite clear from the figure 
|^ that for a practical simulation of the domain- 
wall fermion, one either chooses a large /3 with the 
conventional Wilson action or an improved action 
keeping lattice spacing large. 

To summarize, we have studied the residual 
chiral symmetry breaking present in domain- wall 
fermion by measuring the eigenvalues of the her- 
mitian domain- wall Dirac operator corresponding 
to the topology of the lattice gauge configura- 
tions. Individual eigenvalues for the topological 
zero modes show clear exponential behavior in L s . 
We regard these eigenvalues as the induced mass 
for the surface chiral modes at finite L s separa- 
tion. 

For L s and /?, we see little variation of m e g 
as a function of the volume. This is in some 
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